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ABSTRACT 

Short precession periods like 164d of SS433 can be well determined by obser- 
vations of time scales longer or much longer than the precession period. However, 
it doesn't work for sources with precession periods of millions of years. This paper 
utihzes the particular morphologies of X-shaped sources, so that the 3 dimension 
kinematics of lobes can be obtained. Thus, for the first time, the million-year 
precession period of X-shaped sources by observer on the Earth can be deter- 
mined elegantly: (6.1 ± 1.5)Myr, (1.8 ± 0.5)Myr, and (3.2 ± 1.2)Myr for 3C52, 
3C223.1 and 4C12.03 respectively. The result naturally explains the asymme- 
try displayed in the morphology of these sources, and the effect of propagation 
time on the diversity of morphologies is well demonstrated. The precession pe- 
riod may originate from long-term effects of a binary super-massive black hole 
system, which is a potential source of gravitational wave radiation. 

Subject headings: galaxies: active — galaxies: jets 
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Introduction 



years (ICheung 



Fi rst discovered in 1974 (IHogbom fc Carlsson Ill974l ) and growing rapidly in recent 



20071 ). the peculiar radio morphologies, X-shaped extragalactic radio 
sources, are characterized by two low surface brightness wings oriented at an angle to the 
high surface brightness lobes, giving the total an X-shape. Recently, X-shaped galaxies 
are being considered as potential transition between Fanaroff- Riley (FR) type I and 



II (ILandt et al. 



2010f ). 



Several formation scenarios have been proposed. One is the back-flow of plasm a 



from the active lobes into the w ings fiLeahy fc Williams 



1984 



Hodges-Kluck fc Reynold 



Capetti et al. 



2002 



201 ll ). with subsequent buoyant expansion. It has been argued 



that the expansion of wings is subsonic, and it becomes untenable fo r X-shaped sources 



with wings longer than the active lobes (iDennett-Thorpe et al. 



is the conical precession of the jet axis ( 



Parma et al. 



1985 



2002). T 



le second scenario 



Mack et al. 



19M), which 



implies a scenario of ballistic jet motion plus jet precession, predicting spiral pattern. 
However, this model requires a specific accident of the positions at which the source first 
switched on and its position now. Mo reover, it can not explain the notable asymmetry in 



3C223.1 flDennett-Thorpe et al. 



20021 ). 



Other two explanations have received much attention lately. They are in agreement 
that the wings are the relics left over from a rapid realignment of a central super-massive 
black-hole (SMBH) accretion disk system. The realignment can be a result of a relatiy ely 



recent merger 



Mezcua et al. 



of a s uper-massive binary black hole (SMBBH) ( Merritt fc Ekers 



201 ll ) or due to disk-instability (IDennett-Thorpe et al 



20021 ) 



2002 



Nevertheless, in these two scenarios the morphology of X-sources apparently 
requires a rapid change of jet orientat ion. A number of X-sources have companion 



galaxies ( iDennett-Thorpe et al. 



20021 ) . and the host galaxy of 3C293 shows clearly 
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interaction. Moreover, the double-peaked lo w-ionization emissio n lines in the nucleus of 



20071 ) provide an interesting 



a galaxy associated with X-shaped structure (lZhang et al. 
signature of link between X-source and SMBBH. Such a link makes X-sources potential 
sources of gravitational wave radiation, which leads to a steep surge of interest in them 



recently (jKomossa 



20031 ). although they have been known for decades. 



In fact, both the jet precession and reorientation scenarios agree that the wings are 
relics of the previous jet, and the lobes are produced by jet in action. Thus, regardless of 
the mechanism of the change of the jet axis. We establish a coordinate system, in which 
the 2-dimensional morphology can be fitted by the simplest geometry. And together with 
the constraint imposed by the simultaneous arrival time of photons from the south and 
north lobes of an X-source, the most collimated components of Fig. [H the 3-dimensional 
kinematics of the morphology can be obtained. Consequently, the time scale of the formation 
of the X-shaped morphology (wing and lobe) can be determined i n a general mann er. 



Applying to three X-sources of FR-II, 3C52, 3C 223.1 and 4C12.03 flMerritt fc Ekers 



Dennett-Thorpe et al. 



2002 



Lai fc Rao 



2002; 



20071 ) as shown in Fig. [T], the precession periods 



(by Earth observer) of (6.1 ± 1.5)Myr, (1.8 ± 0.5)Myr, and (3.2 ± 1.2)Myr are obtained 
respectively. Throughout this letter we assume a Qa=0.73, fim=0.27 cosmology with 
/7o=71km s-^Mpc-i. A nd the red-shift, of 3C52, 3C223.1 and 4C12.03 are 0.2 854, 0.1074 



and 0.1570 respectively (lMerritt fc Ekers 



2002 



Saripalli fc Subrahmanyan 



20091). 



2. The New Approach 



The non-ballis tic model (Gong 



motions of AGNs ( Kellermann et al. 



20081) has been used to interpret the non-radi al jet 



2004 



Agudo et al. 



2OO7I : 



Lister et al. 



2OO9I), in 



which a knot can be produced by a continuous jet interacting with ambient matter in 
different directions during the precession of the jet axis. Approximately equal knot-core 
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Fig. 1. — The observed morphology of three X-shaped sources (jLal &: Rao 1120071 ). and the fitted 
traces (ellipses) given by the non-ballistic model. The north is to up, and the east is to left. The 
solid ellipses correspond to the best fit parameters of Table 1, and the dashed ellipses correspond 
to 20% increase in parameters, /, A, and (while holding others constant) for 3C52, 3C223.1 and 
4C12.03 respectively. The red dots represent components of a morphology which the ellipses try 
to fit. The red circles are "radius" of the lobes through which the error in precession phase is 
estimated. 

separation is expected when the power of the jet and matter density of the surrounding 
medium are unchanged in different directions. Such a constant core-knot separation avoids 
the specific accident of the positions required by conical precession model. And due to the 
X-shaped morphologies display similar non-radial characteristic as other AGN sources, it is 
conceivable to apply the non-ballistic model to the X-sources. 

This model can be described by two simple geometric equations. Projecting a knot, i, 
with knot-core separation, R\ to the coordinate system x — y — z, we have. 



i?* [sin A sin / cos rj^ + cos A cos /] 
R" [sin A sin 77*] , 

R' [cos A sin / — sin A cos / cos 17*] 
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where X, I , rj, and R, represent the opening angle of the precession cone, the inchnation 
angle between the jet rotation axis and the line of sight (LOS), the precession phase, and 
the knot-core distance respectively, as shown in Fig. HJ 



LOS-1 



Fig. 2. — Schematic illustration showing an X-shaped source under the non-ballistic precession 
scenario. Different precession states of the X-sources are equivalent to the observation of one source 
at different view angles as LOSi and LOS2, and hence different morphologies. 

In Eq. ([1]) the x-axis is towards the observer. Rotating around the x-axis for angle 
^, so that the new y-axis (A5) will point north, and the new 2;-axis (Aa) will point east. 
Therefore, the coordinate of a feature in the plane of sky is given, 

i?^ = i?* [sin A sin rf] sin ^ + K [cos A sin / — sin A cos / cos rf] cos ^ , 

R\ = i?* [sin A sin ?7*] cos ^ — -R* [cos A sin / — sin A cos / cos ?7*] sin ^ , (2) 

Eq. ([1]) predicts an ellipse, the shape of which is determined by four parameters, ^, A, /, 
and R. Searching in the parameter space of them as shown in the brackets of Table 1, the 
best combination of parameters corresponding to the solid ellipses in Fig. [1] can be found. 

Once the best fit ellipse for a morphology is found, the precession phase of a knot. 
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1]^ = Qt + rjQ, can be given, where Q is the precession velocity of a jet, and rj^ is the initial 
phase of a knot. However, the precession time, t, and the precession velocity, Q, can not 
be separated from r/*. In other words, Q can not be obtained by such a 2-dimensional 
morphology fitting alone. 

Fortunately, a simple constraint can be found to split Q and t, and therefore allows us 
to determine the precession period by 3-dimension kinematic. The time of emission of a 
photon from a knot, temit? can be measured in the reference frame at rest to the core of an 
X-shaped source. This photon can reach an observer on Earth at, tarr = i\mit + djc — xYc, 
where d is the core-observer distance, and x*, which is equivalent to i?^ of Eq. ([T]) is the 
projection of the knot-core separation onto LOS. For simplicity, we can define a time, 

= tarr ~ d/c, SO that the time taken for a photon from a knot to the observer can be 
represented simply by, ta = f^emit — x'^ /c. 

Also measured in the reference frame at rest to the core of an X-shaped source, the 
precession time of a knot can be synchronized to the emitting time, = tg^jj, where is 
given by r/* = VLf + r/g. Thus, the condition of observing the two signals from two opposite 
lobes (2 = 1,2 denote the north and south lobe respectively) at the same time becomes: 



Certain precession phases result in, — /c > ta, which means that the signal has not 
arrived to the observer yet, and is hence unobservable, whereas, f — x^/c < ta means that 
the received signal is from a knot that has been at its emitting site for a period of time. If 
it has afterglow emission that is above the threshold of detection, then it is still observable. 

Therefore, a knot which is unobservable in the case of zero cooling time, f — x^/c < ta, 
becomes detectable provided the emissivity of the knot is above the threshold of detection 
in the cooling time, t*. Such an emission reaches the Earth at: 



f - x'/c = t, 



(3) 



tl + f - X'/C = ta . 



(4) 
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Consequently, with larger and larger separation of the active lobes, a knot in the wings 
corresponds to larger and larger t*, so that it appears more diffused and faint until it 
becomes unobservable. 

Limited by the expansion speeds of low-luminos ity FR-II sources, the advan ce 



20021). Such 



speeds of X-sources are likely not greater than 0.04c (iDennett-Thorpe et al. 
non-relativistic speeds correspond to a negligible Doppler boosting effect, so that the flux 
of a knot depends primarily on the emissivity. Hence, the precession period at the frame 
of the source, Pp, and at the Earth, P°^'^, are related by Pp{l + z) = P°** , where z is the 
red-shift of the source. 

Multiplying the precession velocity, fl, at the two sides of Eq. (|3]), we have, 

rf -rf^- Vtx'/c = ria , (5) 

where rja = fit a- Apparently, an X-source is observed when the photons from the lobes 
arrive at the Earth simultaneously. At this moment the initial phases of the two lobes can 
be treated as, tiq = tiq, without losing generality. Thus ?7q of Eq. ([5]) can be canceled. 

The 77* and of the north and south components can be obtained by fitting the 
2-dimensional X-shaped morphology. The process is simply putting the geometrical 
parameters of Table 1 (except Q and into Eq. (|2]), then find the best-fit parameters 
through minimizing the sum of the square of the residuals of the predicted ellipse from 
the observed morphologies. The pair ellipses (with same group of geometric parameters) 
for each source are required to fit 7-8 components represented by the red dots in Fig. [H 
The solid ellipses of Fig. [1] represent the best fit ones corresponding to parameters shown 
in Table 1, and the dashed ellipses correspond to 20% increase in parameters, /, A, and 
(while holding others constant) for 3C52, 3C223.1 and 4C12.03 respectively. Hence, the 
role of these parameters in the formation of an ellipse is exhibited, i.e., for 4C12.03, the 
20% increase in ^ from its best value changes not only the shape of the dashed ellipse, but 
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also its position, in which case, the morphology cannot be fitted, no matter what other 
parameters are. 

Although making an ellipse through 7-8 points in the morphology strongly constrains 
the fitting parameters, we did find that the morphology of 3C223.1 can be fitted by other 
combination of parameters, i.e., R approximately twice and I half of the corresponding 
parameters of Table 1. However, such a solution predicts much larger discrepancy in 
propagation time between the Northwest (NW) and Southeast (SE) structure than that of 
Table 1, which is contradict to the nearly symmetric structure of 3C223.1, as analyzed in 
Section 3. It is thus excluded. 

Differentiating Eq. ([2]) one has, A_RJ. = ^ fj'^crj, where (Tj{j from 1 to 4), denotes A, 
J, ^ and R respectively, and fj correspond to their partial differentiations respectively. If 
Ai?^, where k represents a and 6 of Eq. ([2]), could be as large as the size of a lobe as shown 
by red circles in Fig. [1] (which is attributed to the error of the precession phase of a lobe), 
then the errors of aj can be obtained by solving four equations, (A_RJ_)^ = J2 fji^'^j)'^^ 
corresponding to two lobes, i = 1,2 (the north and south), of a source. This gives 
conservative errors (up to 43%) to the best fit parameters of Table 1, which considerably 
exceed the 20% parameter errors corresponding to the deviation between the solid and 
dashed ellipses shown in Fig. [H 

Although the 2-dimension morphology fitting perpendicular to the LOS can give 77* and 
x*, the most interesting parameter, Q, can not be extracted from Eq. ([5]). Since the NW and 
SE structure indicates that photons from them must arrive at the Earth simultaneously, this 
can be treated as another constraint (1-dimension), along the LOS. The two active lobes(the 
north and south lobes) with the shortest cooling time always satisfy — x^/c = — x^/c, 
no matter if ^ [x^ x^), or if these values differ largely. Hence the kinematics of the 
two active lobes of an X-source reads, t]^ — Vlx^ /c = rf — Vlx"^ /c, which includes 3-dimension 
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constraints, and from which Q, can be extracted. 

To obtain Q and its error from Eq. both Q and rja can be ordered as variables. 
Thus, Eq. ([5]) corresponds to two hnes {i = 1,2), the cross point of which, rja), represents 
that the two signals arrive on the Earth simultaneously, as shown in Fig [31 With r^* and 
X* of the north and south lobes obtained through morphology fitting (2-dimension), and 
utilize the constraint on photon arrival time (1-dimension), the cross point (fi, can be 
determined by Eq. ([5]), through which the precession velocity, Q, can be constrained into a 
small range elegantly, as shown in Table 1. The error of Q is determined by the errors of 
?7* and x*, where Ax* can be given by the error propagation of parameters in Table 1 via 
Eq. ([1]), and Ar^* is assumed to correspond to the "radius" of a lobe. 



3. Asymmetry 

The fitting parameters of Table 1 indicates that the phase discrepancy of the two 
lobes are ^r/^^ = 185° and 6r]^^ = 189° for 3C52 and 3C223.1 respectively. The ^ 180° 
phase discrepancy corresponds to the observation of an X source through the LOS2 of 
FigO Correspondingly, — x^/c = — x^/c is satisfied in the case x^ ~ x^ and ~ t^. 
As the time and increase, the values of x^ and x^ change similarly. Therefore, the 
emission from the lobes and wings can arrive the Earth at approximately the same time. 
Consequently, the SE and NW structures of these two sources appear comparable in length 
and size. 

In contrast to these two X-sources, the north and south lobes of 4C12.03 correspond to 
a phase discrepancy of dr]^'^ = rj^ ~ rf —270° instead of ~ 180°, by the fitting parameters 
of Table 1. This means that the photons from the active south lobe, which should differ 
by approximately ~ 180° to the precession phase of the north lobe, have not reached the 
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Fig. 3. — The determination of the precession period of three X-sources. The two Unes that cross 
at {n,r]a), represent the same time of arrival of photons from the north and south lobes. The 
dashed lines correspond to uncertainty in r/* and x*, estimated by the size of the lobes. The jet 
precession velocity of a source, il, is thus constrained in the small regions, ABCD. 
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observer on the Earth yet (t^ — x^/c > ta), ahhough such a lobe is observable at the core 
of this X-source. This corresponds to the observation of such a source through the LOSi of 
Fig[2J Therefore, —x^/c = t"^ —x^jc is satisfied in the case a;^ and ii} and also) differ 
significantly. Such a discrepancy in x^ and x^ results in a Southwest (SW) pattern compared 
to that of Northeast (NE). Hence, the significant asymmetry of 4C12.03 is well understood. 

The north and south lobes of 4C12.03 obviously differ in size. By the fitting parameters 
of Table 1, the active north lobe posses the maximum and minimum 77^ (due to i7 < 
in 4C12.03), which means that knots of t > and rj < t]^ don't exist in the NE pattern at 
all. This is in agreement with both observers on the Earth and at the core of this X-source. 
Consequently, the emission of this lobe region can only be extended by the past (cooled) 
components, with t <t^ and t] > rj^. 

Contrarily, for the SW pattern, its emission can be extended by both the past 
components with — x^/c < ta, and by some "future" components with t"^ — x'^/c > ta- 
Because components with t > and r] < iff do exist near the "south lobe" (which would 
be observable at the core of this X-source), and the cooling emission of such "future" knots 
can contribute to the brightness of the "south lobe" as well hy + — x'^j / c = ta (where x / 
denotes the "future" emission site). Therefore, the one way extension of the NE lobe and 
the two way extension of the "SW lobe" lead to a larger "south lobe" than that of the NE 
one. 

To analyze the fine asymmetry in these X-sources, Fig |4] is introduced, which is 
obtained as follows. Putting the obtained r/* and through the fitting of morphologies of 
three sources of Fig [1] into Eq. ([5]), then ordering VL = const, the phase corresponding to 
the time of arrival, rja, versus the phase of precession, 77, can be obtained as shown in FigHJ 
which is actually the evolution of ta and t of the two active lobes. 

As shown in FigHl each source has its own rj — r]a curve and horizontal dashed line, rja. 
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Fig. 4. — The precession phase versus the phase of arrival time, with = const. The cross points 
of the horizontal line with two curves, rj^ and rj^ corresponds to the precession phases of the north 
and south lobe respectively. The cross points of the curves with the vertical lines (r/*') correspond 
to the precession phase of the tail of wings. 
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which cross with a misahgnment angle. In both 3C52 and 3C223.1, the misahgnment angles 
of the south lobe are smaller than those of the north lobe. Moreover, the discrepancy in the 
misalignment angles of the north and south lobe is more obvious in 3C223.1 than in 3C 52. 
This explains the deviation in surface brightness of the NW and SE lobe regions in 3C223.1, 
as shown in Fig [H Because the region near the south lobe is closer to the horizontal line, 
rja, than that of the north one, which corresponds to a shorter cooling time. Hence the 
linking of the south lobe and its neighboring wing region has higher surface brightness than 
those of the north one. 



4. Discussion 



By the fitting of Fig ([T]) and parameters of Table 1, precessing across the lobe- wing 
region takes 1.5 Myr, 0.8 Myr and 1.4 Myr for 3C52, 3C223.1 and 4C12.03 (NE pattern) 
respectively, the time scale of which corresponds to the cooling time discrepancy between 
the active lobe and the tail of wing of these sources. Interestingly, the time scale is 
consistent with the firm upper limits on the particle ages of 34 Myr for 3C223.1, and the 
estimation of the time scale of reorientation of j et axis of no more than a few Myr based on 



spectral gradient ( jPennett-Thorpe et al. 



20021). 



Beside the temporary processes such as merger of SMBBH or disk instabil- 



ity ( iMerritt &: Ekers 



2002 



Dennett-Thorpe et al. 



20021 ) . the change of j et orientation 



disp 



ayed in X-shaped morphology can also be originated from binary effect (IBegelman et al. 



19801 ). which can be either relativistic geodetic precession or Newtonian-driven jet pre- 



cession (jKatz 



jital period of 20 ye ars, and with a 



20111 ) and companion 



19971 ). A SMBBH system with orl 
typical X-source black hole mass of 5 X IO^Mq flMezcua et al. 
mass of 1 X IO^Mq, predicts a precession period of 0.6Myr by the geodetic effect; and a 
precession period of l.OMyr (with disk radius of 10 schwarzschild radius) by the Newtonian 
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driven effect. Consequently, the ~Myr time scale displayed in the three X-sources can 
be well interpreted by either of the two binary effects. This provides another evidence 



to the link between X-sources and SMBBHs. M oreover, the correlat i on of black 



lole 



mass with X-ray luminosity and radio luminosity ( iMerloni et al. 



2003 



Falcke et al. 



20041 ): 



as well as characteristic time-scale of the X-ray variability finding in AGNs and X-ray 



binaries ( iMcHardy et al. 



2006 



Kording et al. 



2007 



Mirabel 



20061 ) suggests that black 



hole physics likely scales with mass. The ~Myr precession periods of the three X-sources 
revealed here further support such a correlation with black hole mass. 

In the non-ballistic precession model, the diversity of morphology of X-sources can be 
simply understood by the received photons from emission sites at different precessi on cones 



The t echnique is applicable to other X-sources, especially the X-shaped candidates ([Cheung. 
20071 ) in the future. 
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Table 1: The parameters extracted by morphology fitting and 3-dimension kine- 
matics. 

Source C(0, 27r) 7]^(0,27r) ry2(0,27r) /(0,7r/2) Ai(0,7r/2) i?(0, 2) n(arcsec/yr) n°'"'(arcsec/yr) 
3C52 5.2 ±0.14 0.88 ± 0.44 -2.4 ± 0.34 0.52 ± 0.065 0.43 ±0.040 0.40 ± 0.093 0.27 ±0.062 0.21 ±0.048 
3C223.1 5.1 ±0.061 1.1 ±0.50 -2.2 ± 0.41 0.62 ± 0.23 0.45 ± 0.077 0.15 ± 0.021 0.83 ± 0.23 0.75 ± 0.21 
4C12.03 3.9 ± 0.046 -1.5 ±0.66 3.1 ± 0.50 0.65 ±0.15 0.35 ± 0.079 0.35 ±0.11 -0.47 ±0.17 -0.41 ±0.15 

The precession velocities, 0,, in the frame of the source, are inferred from Fig O which is related 
to the frame of the Earth hy Q = Q°''^{1 + z). The rest are precession parameters obtained by the 
fitting of Fig[TJ Angular parameters are in rad, in which r/^ and ry^ represent the phase of the north 
and south lobe respectively. The opening angle of precession cone, A^, corresponds to the NW(NE 
for 4C12.03) pattern, and the opposite pattern is fitted by = vr — A^. The distance R is in Mpc, 
which is obtained by multiplying the arcsec value (through morphology fitting) by the angular size 
distance (through the cosmology parameters and the measured redshift). 



